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CROSS REFERENCES TO RELATED APPLICATIONS 

The present application is a continuation-in-part of co-pending U.S. 
Patent Application Serial No. 09/766,879, filed January 19, 2001. The entire 
contents of the above application is incorporated herein by reference. 

5 BACKGROUND OF THE INVENTION 

Adenocarcinoma of the lower esophagus develops almost exclusively in patients 
with Barrett's esophagus (BE), a condition characterized by the presence of metaplastic 
columnar epithelium. While the prognosis of patients diagnosed with adenocarcinoma 
is poor, the chances of successful treatment increase significantly if the disease is 
10 detected at the dysplastic stage. The surveillance of patients with BE for dysplasia is 
challenging in two respects. First, dysplasia is not visible during routine endoscopy. 
Thus, numerous random biopsies are required. Second, the histopathologic diagnosis of 
dysplasia is problematic, as there is poor inter-observer agreement on the classification 
of a particular specimen, even among expert gastrointestinal pathologists. Optical 
15 techniques, such as fluorescence, may significantly enhance the endoscopist's ability to 
detect these early dysplastic changes in BE. Indeed, fluorescence spectroscopy studies 
using exogenous fluorophores, such as Photofiin® and aminolevulinic-acid induced 
protoporphyrin IX, show that there is a significant difference between the measured red 
fluorescence of the carcinomatous and non-dysplastic tissue as a result of the 
20 preferential accumulation of the drug. Initial autofluorescence spectroscopy studies 
performed at 410 nm excitation report promising results for detecting high-grade 
dysplasia. However, focal high-grade and low-grade lesions could not be detected 
reliably. Thus a continuing need exists for further improvements in the optical 
measurements used to detect early stage carcinomas. 

25 SUMMARY OF THE INVENTION 

The present invention relates to a combination of spectroscopic systems that can 
improve the sensitivity and accuracy of measuring the state of tissue. For example, a 
plurality of different spectroscopic systems can be used to determine whether a region of 
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tissue is normal, pre-cancerous or cancerous. In a preferred embodiment dysplasia 
detection in patients with Barrett's esophagus (BE) can be performed. In a preferred 
embodiment, fluorescence, reflectance and light scattering spectroscopies provide 
complementary information about the biochemical, architectural and morphological 
5 state of tissue and the corresponding changes that occur during the progression of 

disease, and in particular, of dysplasia. In a preferred embodiment, the system has been 
developed for combining these three methods to provide for detection, mapping, and/or 
imaging of tissue. In particular, methods of processing the spectral or image data with a 
data processor programmed to perform an automated diagnostic routine is useful to 
1 0 provide the results of a plurality of different measurement modalities including 

fluorescence and light scattering spectroscopy. A further preferred embodiment of the 
invention utilizes a tri-modal system to guide a biopsy procedure. 

Of importance in this system for real-time measurements is the simultaneous or 
near simultaneous collection of light from the same spot or region of interest. For 
1 5 example, the collection of light includes the collection of a measured fluorescence and 
reflectance spectra. The detected reflectance spectrum is processed to remove a 
diffusive background component. Once this is accomplished, then a component of the 
light that is periodic in wavelength is measured. This component arises from the light 
that is Mie-scattered by surface epithelial cell nuclei, for example. By analyzing the 
20 amplitude and frequency of the periodic structure, the density and size distribution of 
these nuclei can be extracted. For the reflectance, light scattered and fluorescence 
components to be properly correlated and used to assess a given region of interest, there 
is preferably substantial overlap of the excitation light for both the reflectance and 
fluorescence measurements. The reflected light is used to both correct the fluorescence 
25 spectrum and to generate a light scattered spectrum based on the use of the periodic 
structure contained therein. The apparatus delivers both excitation components to the 
region of interest through the same distal surface of the probe, preferably through the 
same optical fiber or collection of fibers. 
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The biopsy channel of an endoscope can be used to insert the fiber optic light 
delivery and collection system used to obtain measurements. Alternatively, a small 
diameter endoscope, 5 mm or less in diameter, for example, can include the light 
delivery and collection system suitable for many applications. A preferred embodiment 
5 of the system can use a single fiber system for delivery and collection, or alternatively, 
can employ a central delivery fiber and six collection fibers concentrically arranged 
around the delivery fiber. The proximal end of the light delivery and collection probe is 
optically coupled to both a broadband flash lamp and a monochromatic source such as a 
laser. A rotating filter or dye wheel can be used to rapidly switch the excitation 
CI 1 0 wavelength over a selected range. 

U! The need for using reflected light arises from the need to correct for the effects 

| S of hemoglobin absorption on the measured integrated tissue fluorescence intensity. The 

B combination of fluorescence and reflectance spectroscopies can be applied to remove 

distortions introduced by scattering and absorption into the entire measured tissue 
2 1 5 fluorescence spectrum. The undistorted or modified fluorescence spectrum can serve as 

a sensitive indicator of tissue biochemistry, while reflectance and light scattering 
O representations provide morphological information on tissue architecture and epithelial 

cell nuclei. The present invention can include the simultaneous use of all three 
spectroscopic methods for characterizing tissue and diagnosing disease. The following 
20 demonstrates that the combined use of all three techniques provides improved results as 
compared to the results of each technique individually, in terms of detecting not only 
high-grade, but also low-grade dysplastic changes in BE, for example. 

In a preferred embodiment, a method of analyzing spectral data to measure a 
structure in a layer of tissue includes providing a light collection system to collect 
25 fluorescent and reflected light from the tissue at a plurality of wavelengths and detect 
the collected light. The method further includes forming a fluorescence representation 
and a scattered light representation as a function of wavelength from the detected light. 
The method includes determining a characteristic of the tissue layer with the fluorescent 
representation and the scattered light representation. The method includes using the 
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scattered light representation to determine a size of a structure, for example, the nuclei 
of cells, within the tissue layer. A diffuse reflectance spectrum can be analyzed in a 
preferred embodiment. The method, in the preferred embodiment includes generating a 
lookup table with different sizes of scatterers. The fluorescence representation can be 
5 corrected to obtain an intrinsic fluorescence spectrum. 

In a preferred embodiment of the present invention a method for analyzing tissue 
spectra includes acquiring a fluorescence spectrum, acquiring a reflectance spectrum 
and then processing the two to provide an intrinsic fluorescence spectrum, a diffuse 
M: reflectance spectrum and a light-scattering spectrum. The method includes determining 

S 10a biophysical tissue characteristic from the intrinsic fluorescence spectrum, the diffuse 

!" reflectance spectrum and the light scattering spectrum. The method may include 

U 1 performing a discrimination analysis using the diagnostic parameters from the diffuse 

reflectance spectrum, the diagnostic parameters from the intrinsic fluorescence spectrum 
and the diagnostic parameters from the light scattered spectrum and then providing a 

L„.! 

15 diagnosis. 



l"""r. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 schematically illustrates a preferred embodiment of a system for 
performing measurements in accordance with the invention. 

Figure 2A illustrates graphically fluorescence from a non-dysplastic Barrett's 
20 esophagus site, 337 nm excitation wherein the measured spectrum is the solid line; 
extracted intrinsic fluorescence is the dashed line. 

Figure 2B illustrates graphically the corresponding reflectance spectrum in 
accordance with the fluorescence spectrum illustrated in Figure 2A. 

Figures 3A-3D illustrate the fluorescence spectra from non-dysplastic (solid 
25 lines), low-grade (dashed-lines) and high-grade dysplastic (dotted lines) Barrett's 

esophagus (BE) sites. The measured and corresponding extracted intrinsic fluorescence 
for excitation at 337 nm ((A) and (B)) and 397 nm ((C and (D)) are shown. Spectra are 
normalized to their peak intensities. Note the significant lineshape changes. The mean 
± standard deviation values are displayed for each category. 
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Figure 4 illustrates the scores of three principal components extracted from the 
decomposition of intrinsic fluorescence spectra at 337, 397 and 412 nm excitation used 
to distinguish high-grade dysplasia (diamonds) from non-dysplastic and low-grade 
dysplasia (squares) Barrett's esophagus (BE) sites. At 337 nm excitation, 
5 decomposition was performed in the 460 to 520 nm region of the intrinsic fluorescence 
spectra, as this is the wavelength range within which spectral differences are most 
pronounced. Similarly, at 397 nm excitation principal components were extracted from 
the intrinsic fluorescence spectra between 600 and 650 nm. PC 1 indicates the first 
,jl principal component and PC 3 the third. 

y 1 0 Figure 5 illustrates the scores of three principal components extracted from the 

H decomposition of the entire intrinsic fluorescence spectra at 337, 397 and 412 nm 

jl excitation used to distinguish dysplastic (low- and high-grade; diamonds) from non- 

dysplastic (squares) Barrett's esophagus (BE) sites. PC 2 indicates the second principal 
I component. 

;ssy 

15 Figure 6 is a reflectance spectrum of a non-dysplastic Barrett's esophagus site. 

I The solid line represents the measured data and the dashed line represents the projected 

• ; 3 properties of the tissue based on known scattering and absorption properties of the 

tissue. 

Figure 7 illustrates graphically the slopes and intercepts of wavelength 
20 dependent tissue reduced scattering coefficient, ju s ' (X) , for non-dysplastic (squares), 
low-grade (diamonds) and high-grade (triangles) dysplastic Barrett's esophagus (BE) 
sites. In each case, a straight line was fit to ju s 1 (X) , extracted from the reflectance 
spectrum. A log-log scale is used to facilitate visualization of all the data points. 

Figure 8 illustrates the total number of nuclei per mm 2 plotted as a function of 
25 percentage of enlarged nuclei (diameter > 10 ^m), as determined from the 

light-scattering model analysis. Non-dysplastic Barrett's are represented by squares, 
low-grade dysplasia are represented by circles; and high-grade dysplasia are represented 
by triangles. 
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Figures 9A-9D illustrate graphically the measured and intrinsic fluorescence 
spectra of normal squamous epithelium (solid line), benign biopsied sites (dashed lines) 
and high-grade SILs (dotted lines). Curves A and B show differences in the lineshape 
and intensity observed at 337 nm excitation. Curves C and D show intensity differences 
5 observed in the intrinsic fluorescence excited at 358 nm. 

Figure 10 illustrates relative NADH fluorescence levels as a function of relative 
collagen fluorescence levels plotted for normal squamous epithelium (squares), benign 
biopsies (circles) and HSILs (triangles). The three arrows point to the benign biopsies 
that were classified as "mature squamous epithelium." The remaining benign biopsies 
10 were classified as "squamous metaplasia." 

Figure 1 1 A shows reduced scattering coefficient at 400 nm for benign biopsied 
and high-grade dysplatic sites. 

Figure 1 IB shows nuclear crowding or total number of nuclei per mm 2 plotted as 
a function of the standard derivation of the nuclear size population for a particular site. 
1 5 The line is determined by logistic regression. 

Figure 12 illustrates a process sequence in accordance with a preferred 
embodiment of the present invention. 

Figure 13 illustrates a flow chart of a preferred embodiment for measuring a 
structure in a layer of tissue in accordance with the present invention. 
20 Figure 14 illustrates the details of the step of calibration referred to in Figure 13 

in accordance with a preferred embodiment of the present invention. 

Figure 15 is a flow chart of a preferred embodiment of the present invention 
detailing the analysis of the diffuse reflectance spectrum (DRS) in accordance with the 
present invention. 

25 Figure 16 is a flow chart detailing the analysis of the intrinsic fluorescence 

spectrum (IFS) in accordance with a preferred embodiment of the present invention. 

Figures 17A and 17B illustrate a flow chart detailing the analysis of the light 
scattering spectrum (LSS) in accordance with a preferred embodiment of the present 
invention. 
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Figure 18 illustrates the tri-modal system analysis in accordance with a preferred 
embodiment of the present invention. 

Figure 19 illustrates an alternate preferred embodiment of a tri-modal system 
analysis in accordance with the present invention. 
5 The foregoing and other objects, features and advantages of the invention will be 

apparent from the following more particular description of preferred embodiments of 
the invention, as illustrated in the accompanying drawings in which like reference 
characters refer to the same parts throughout the different views. The drawings are not 
H= necessarily to scale, emphasis instead being placed upon illustrating the principles of the 

p 10 invention. 

j }j DETAILED DESCRIPTION OF THE INVENTION 

— Measurements are performed, for example, using a fast excitation-emission 

O matrix (EEM) instrument 10 that has been described in greater detail in U.S. 

Li. 

L Application No. 09/238,664, filed January 26, 1999, the entire contents of which is 

55 15 incorporated herein by reference. The excitation light source of this fast-EEM system 
III can include a nitrogen laser 12 emitting at 337 nm (Laser Science, Inc., Franklin, MA; 

Model: VSL-337MD) pumping 10 dye cuvettes precisely mounted on a rapidly rotating 
wheel 16. In this manner, in a preferred embodiment, eleven different excitation 
wavelengths are obtained between 337 and 620 nm and coupled using optical system 18 
20 into the delivery fiber of a 1 mm diameter optical fiber probe 20. For the reflectance 
measurements, white light (350-700 nm) from a Xe flash lamp 14 (Perkin Elmer 
Optoelectronics, Salem, MA) is coupled into the same probe. Alternatively, for other 
embodiments involving measurements in the bladder, for example, a XeCl excimer laser 
emitting at 308 nm can be used. In a preferred embodiment, the probe is composed of 
25 six collection fibers 28 surrounding the central light delivery fiber 26, and it is covered 
with a protective, transparent optical shield at the distal end 22 as shown in Figure 1. 

During endoscopy, the probe is inserted into the accessory channel of the 
endoscope and brought into gentle contact with the tissue, thus providing a fixed 
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delivery-collection geometry. The reflected and fluorescence light is collected by the 
probe and coupled to a spectrograph 30 and detector 22. A second synchronized wheel 
36 is used to block the fluorescence excitation wavelength. The average of three sets of 
spectra from each site is used for analysis using a data processing system 34. 

5 Immediately following data acquisition, the probe is removed and a slight demarcation 
remains on the tissue for 30 to 60 seconds as a result of the probe contact. This 
endoscopically apparent marker is used as a guide for taking a biopsy from the same site 
at which spectra were acquired. The biopsy specimen is interpreted and classified by a 
gastrointestinal pathologist. If a dysplastic lesion is suspected, the specimen is reviewed 

10 and the diagnosis confirmed by a second gastrointestinal pathologist, in accordance with 
the standard of care. Data has been analyzed from 26 non-dysplastic Barrett's esophagus 
sites (9 patients), 7 low-grade (4 patients) and 7 high-grade (5 patients) dysplastic sites. 

Three types of spectroscopic information are acquired, preferably in less than 
one second in a preferred embodiment. Fluorescence spectra at eleven different 

15 excitation wavelengths, reflectance spectra and light scattering spectra are obtained. 
Each type of spectrum can be analyzed in a manner that provides information about 
biochemical and morphological changes that occur during dysplastic transformation. 
Fluorescence spectroscopy can provide valuable information about changes that take 
place in tissue biochemistry during the development of dysplasia. However, the 

20 measured tissue fluorescence spectra can be distorted significantly by unrelated 

scattering and absorption events. To remove these distortions, the fluorescence spectra 
are analyzed in combination with information from the corresponding reflectance 
spectra. The success of this procedure is predicated on the fact that the fluorescence and 
reflectance spectra collected from a specific site or region of interest using the same 

25 light delivery/collection geometry undergo similar distortions. By extracting the 
intrinsic (undistorted) tissue fluorescence, changes in tissue biochemistry may be 
isolated in a more sensitive and specific manner. 

Principal component analysis and logistic regression is employed to determine 
the correlation between spectral features of the intrinsic fluorescence and 
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histopathological diagnosis. To analyze this relatively small data set in an unbiased 
manner, "leave-one-out" cross-validation can be used. Specifically, the principal 
components of the intrinsic fluorescence spectra that described the spectral features that 
change during the progression of dysplasia are selected. The corresponding scores (the 
5 coefficients describing the contributions of the principal components to the overall 
spectra) can be used to determine the ability of the system to distinguish (a) high-grade 
dysplasia from low-grade dysplastic and non-dysplastic Barrett's esophagus (BE), and 
(b) dysplastic (low- and high-grade) from non-dysplastic Barrett's esophagus (BE). To 
f achieve that in an unbiased manner, the following procedure can be performed. The 

□ 10 scores from a particular site are eliminated, and logistic regression is used to form a 
i I decision surface that classifies the remaining sites in a manner that optimizes agreement 

with the histopathological classification. The resulting decision surface maybe then 
used to classify the excluded site. This process can be repeated for each of the sites. 
This method provides use of a relatively small data set to validate the performance of a 
1 5 decision surface without bias. The decision surface varied minimally during this 

procedure, indicating the reliability of the technique. Sensitivity and specificity values 
can be determined by comparing the spectroscopic classification with histopathology. 
Statistical analysis maybe performed using Matlab statistic software (The Math Works, 
Inc, Natick, MA) in accordance with a preferred embodiment. 
20 The measured reflectance spectra can be analyzed using a representation based 

on diffusion theory which expresses the reflected light as a function of the absorption 
(ju a (X)) and reduced scattering (ju s '(X)) coefficients of tissue. This analysis provides 

information about the architecture and morphology of mainly the connective tissue, i.e. 
the lamina propria and the submucosa, as the collected light originates within 500-700 
25 ]im from the tissue surface. The diagnostic value of the resulting tissue scattering 

coefficient values can be determined by correlating the results of logistic regression and 
cross-validation with histopathological classification, as in the intrinsic fluorescence 
case. 
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A small fraction (2-5%) of the detected reflected light originates from light 
collected by the probe after being scattered only once by the tissue. This method 
described generally herein as light scattering spectroscopy is described in greater detail 
in U.S. Patent No. 6,091,984, issued on July 18, 2000, the entire contents of which is 

5 incorporated herein by reference. Additional methods for measuring tissue structure are 
described in International Application No. PCT/US98/2145, filed on October 9, 1998 
designating the United States, the entire contents of which is also incorporated herein by 
reference. The intensity of this singly-scattered light contains a component which is 
periodic in inverse wavelength, the magnitude and frequency of which depends on the 

10 number and size of the nuclei in the epithelial cell layer. This periodic signal is 

analyzed to determine the number and size of epithelial cell nuclei corresponding to a 
particular site. Logistic regression and cross-validation are then used to compare the 
spectroscopic classification with that of histopathology. To optimize sensitivity and 
specificity, the posterior probability threshold for separating high-grade dysplasia from 

15 non-high-grade dysplasia sites is set to 0.3 in one preferred embodiment. 

Finally, results from all three spectroscopic techniques are combined to 
determine whether the number of correctly classified sites can be improved. 
Specifically, a site is assigned a classification that is consistent with results from at least 
two of the three analysis methods, and this classification is compared to histopathology. 

20 Figure 2 A shows a typical fluorescence spectrum excited with 337 nm light from 

a non-dysplastic Barrett's esophagus (BE) site (solid line). There are two peaks, which 
can be attributed to the presence of two different tissue fluorophores. Note that the 
fluorescence intensity decrease between these two peaks occurs at the wavelength range 
in which hemoglobin absorbs light very efficiently. The effects of hemoglobin 

25 absorption are observed in the corresponding reflectance spectrum, which exhibits 
minima at approximately 420, 540 and 580 nm, corresponding to oxy-hemoglobin 
absorption peaks (Figure 2B). When the measured fluorescence spectrum of Figure 2 A 
is processed in combination with the corresponding reflectance spectrum of Figure 2B 
as discussed herein, the intrinsic (undistorted) tissue fluorescence spectrum at the 
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particular excitation wavelength is obtained (Figure 2A ? dashed line). Note that this 
spectrum consists of a single broad peak. 

Significant differences are observed in the representation of intrinsic tissue 
fluorescence of non-dysplastic and dysplastic Barrett's esophagus (BE) sites excited at 
5 337 nm (Figures 3A and 3B) and 397 nm (Figures 3C, 3D). At 337 nm excitation the 
lineshape of the dysplastic sites broadens and shifts to the red region of the spectrum 
during the progression from non-dysplastic, to low-grade, to high-grade dysplasia. At 
397 nm excitation, the fluorescence increases in the red region of the spectrum for the 
u dysplastic BE sites. Similar changes are observed at 412 nm excitation. 

1 0 These differences can be used to develop algorithms for detecting dysplasia in 

□ 

HI BE. Specifically, principal component analysis, logistic regression and leave-one-out 

j -7" cross-validation are employed to determine the sensitivity and specificity with which 

W one can separate (a) non-dysplastic from dysplastic (low- and high-grade) tissue, and (b) 

high-grade dysplasia from low-grade and non-dysplastic BE epithelium. In each case, 
y- 1 5 the scores of one of the first three principal components extracted from the intrinsic 
i! fluorescence spectra at 337, 397, and 412 nm excitation are used (Figures 4 and 5). The 

Q selected principal components describe the observed spectral differences. From this 

fll 

analysis, sites with high-grade dysplasia can be differentiated from low-grade and 
non-dysplastic sites with high sensitivity and specificity (Table 1). Additionally, 
20 dysplastic and non-dysplastic epithelia can be distinguished with very high sensitivity 
and specificity. 

Table 1 illustrates the accuracy of spectroscopic classification of non-dysplastic 
(NDB), low-grade (LGD) and high-grade dysplastic (HGD) tissue in Barrett's 
esophagus. 
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Table 1 





HGD vs. (LGD and 
NDB) 


(LGD and HGD) vs NDB 




Sensitivity 


Specificity 


Sensitivity 


Specificity 


Intrinsic fluorescence (IF) 


100% 


98% 


71 % 


92% 


Reflectance (R) 


86% 


100% 


79% 


88% 


Light Scattering (LS) 


100% 


91% 


93% 


96% 


Combination of IF, R and LS 


100% 


100% 


93% 


100% 



Reflectance spectra can be analyzed using a mathematical representation to 
obtain detailed information about the scattering and absorption properties of the bulk 
H tissue. A typical reflectance spectrum with the corresponding fit obtained using this 

g representation is shown in Figure 6. 

% ? 5 This type of analysis shows that the reduced scattering coefficient ju s ' (X) of 

Barrett's esophagus tissue changes gradually during the progression from non-dysplastic, 
W to low-grade, to high-grade dysplasia. For example, at 400 nm the ju s f (X) of high-grade 

O dysplastic tissue (1.3 ± 0.2 mm" 1 ) is lower than that of low-grade dysplastic tissue (1.8 ± 

0.3 mm" 1 ) which, in turn, is lower than that of non-dysplastic Barrett's esophagus (BE) 
S 10 tissue (3 ± 1.6 mm" 1 ). Additionally, the wavelength dependence of ju s ' (X) changes 
■ ' during the development of dysplasia. To describe these changes, a straight line is fit to 

ju $ ' (X) and the intercept and slope of thus line are used as diagnostic parameters (Figure 
7). Using logistic regression and leave-one-out cross-validation, the sensitivity and 
specificity for 10 classifying tissue in accordance with histopathology are determined. 
15 This method results in slightly lower overall sensitivity and specificity values than those 
achieved with the intrinsic fluorescence spectra (Table 1). 

The reflectance spectra are further processed in a manner that allows extraction 
and analysis of the backscattered light from the epithelial cell nuclei. The results of this 
analysis are displayed in Figure 8. The ordinate of Figure 8 represents the number of 
20 nuclei per square mm, indicative of the degree of nuclear crowding, and the abscissa 
represents the percentage of enlarged nuclei, defined as nuclei having diameter greater 
than 10 microns. Note that the non-dysplastic samples are concentrated in the lower 
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left-hand corner, indicating cell nuclei that are small and free of crowding. As dysplasia 
progresses, the data points move to the upper right, indicating nuclear enlargement and 
crowding, in agreement with the findings of histopathology. This technique is superior 
in terms of separating the dysplastic (low- and high-grade) from the non-dysplastic BE 

5 sites (Table 1). 

The ability to characterize dysplastic and non-dysplastic tissue in BE is 
improved by combining the information provided by each one of the spectroscopic 
techniques, obtained simultaneously with our system. When a spectroscopic 
classification is consistent with at least two of the three analysis methods, high-grade 

1 0 dysplasia is identified with perfect sensitivity and specificity, and dysplastic tissue is 

distinguished from non-dysplastic tissue with perfect specificity, while maintaining very 
high sensitivity (Table 1). 

Spectroscopic techniques use information contained in light signals to assess the 
state of biological tissue. Optical fiber technology allows spectroscopic methods to be 

15 applied as a diagnostic tool for a wide range of tissues that are accessible 

endoscopically. The use of spectroscopy as described herein can be used for improving 
the physician's ability to detect pre-cancerous (dysplastic) and early cancerous lesions in 
many organs, such as the oral cavity, the cervix, the lung, the breast and the 
gastrointestinal tract. Depending on the technique employed, specific information can 

20 be acquired about tissue biochemical, architectural and morphological features. 

Microscopic changes in these features that occur during the progression of dysplasia 
may be detectable spectroscopically before the manifestation of macroscopic changes 
that are visible endoscopically. Additionally, spectroscopic techniques are 
non-invasive, allowing study of the tissue in its native state, free of artifacts introduced 

25 by cutting and chemically processing the tissue. In principle, spectroscopic signals can 
be analyzed in real-time, thus guiding the physician to biopsy areas that are likely to 
yield significant pathology or possibly allowing her to make an immediate decision on 
the type of intervention that is required for successful treatment of the patient. 
Furthermore, the spectroscopic signals carried by light can be used as objective guides 
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for assessing a particular tissue site, especially in areas in which the intra- and 
inter-observer agreement on the classification of disease is not very good. 

The targets of fluorescence spectroscopy include tissue biochemicals such as 
NADH, FAD, collagen, elastin and porphyrins. Exogenous or exogenously-induced 
5 chromophores that have been shown to accumulate preferentially in the diseased areas 
can also be used. The detection of high-grade dysplasia using tissue autofluorescence 
excited at 410 nm have been conducted. The difference between the measured 
integrated intensity-normalized fluorescence and the mean normalized fluorescence 
5f from normal esophageal tissue was used for the diagnostic procedure. The main 

u 1 0 spectral features that resulted in good differentiation between high-grade dysplastic and 
S j non-dysplastic tissues were the presence of decreased fluorescence around 470-480 nm 

: and increased fluorescence in the red region of the spectrum for the high-grade 

Ly dysplastic tissues. However, this process does not classify correctly sites with low- 

□ grade or focal high-grade dysplasia. 

[7 15 In the present example illustrating an embodiment of the invention fluorescence 

33 spectra at 1 1 different excitation wavelengths between 337 and 610 nm were obtained. 

f|j Thus, instead of a single fluorescence spectrum an excitation-emission matrix (EEM) is 

collected. EEMs can be used to identify the excitation wavelengths at which tissue 
classification is optimized. Additionally, EEMs can assist in identifying the origins of 
20 the measured fluorescence signals in a more reliable manner. Nevertheless, as shown in 
Figures 2 and 3, these measured EEMs can be distorted significantly by tissue scattering 
and absorption. To eliminate artifacts introduced by changes in scattering or absorption, 
rather than by tissue biochemistry, corresponding reflectance spectra can be used which 
are affected in a similar manner by scattering and absorption events. Once the distorted 
25 measured tissue fluorescence spectra are rectified using the reflectance, tissue 

fluorescence excited at 337 nm broadens and shifts to longer wavelengths in a very 
consistent manner as the tissue progresses from non-dysplastic to low-grade to 
high-grade dysplasia (Figure 3). These spectral changes are consistent with the presence 
of increased NADH levels in dysplastic tissue. Our findings at 397 and 412 nm 
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excitation are attributed to endogenous porphyrins. The spectra corresponding to the 
high-grade dysplasia sites appear slightly distorted around 470 nm, even after correcting 
for the effects of scattering and absorption. This suggests that this difference arises as a 
result of biochemical changes rather than absorption changes. 
5 To demonstrate the level of significant changes that are observed in tissue 

fluorescence during the development of dysplasia, the scores of one of the first three 
principal components which described over 99% of the variance observed in the 
intrinsic fluorescence spectra excited at 337, 397 and 412 nm can be used. 
Subsequently, one can use logistic regression and leave-one-out cross-validation to 
1 0 estimate and validate in an unbiased manner the sensitivity and specificity with which 
one can distinguish (a) high-grade dysplasia from low-grade and non-dysplastic tissue, 
and (b) dysplastic (low- and high-grade) from non-dysplastic tissue to separate 
high-grade dysplasia from low-grade and non-dysplastic tissue, spectroscopic 
classification is consistent with histopathology in all but one case. Additionally, one 
1 5 can distinguish dysplastic from non-dysplastic tissue with very high sensitivity and 
specificity. 

Reflectance spectroscopy can be used not only to remove the distortions 
observed in the measured tissue fluorescence spectra, but also to provide very detailed 
and potentially useful information about morphological and architectural features of the 
20 tissue. Specifically, as shown in Figure 6 the observed tissue reflectance spectra can be 
used in terms of two parameters that are determined by tissue scattering and absorption. 
For example, changes in the concentration or the oxygen saturation of hemoglobin, the 
main absorber in the visible spectrum for this tissue type, result in concomitant changes 
in the absorption coefficient of tissue. Alterations in the architecture of the connective 
25 tissue collagen fibers, one of the main contributors of tissue scattering, lead to a 
modified tissue scattering coefficient. The analysis suggests that the scattering 
coefficient of tissue decreases significantly during the development of dysplasia, 
suggesting that changes that are not observed histopathologically are taking place within 
the lamina propria and submucosa before the onset of invasion. Recently, it has been 
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shown that an increased level of cysteine and serine proteases is found in gastric and 
colorectal cancerous and precancerous lesions. The findings related to the decrease in 
the value of the scattering coefficient during the progression of dysplasia are consistent 
with the presence of such enzymes, which can result in a less dense collagen matrix, for 
5 instance. The change in the slope of ju s ' (A) as a function of wavelength suggests that 
the mean size of the tissue scattering particles is changing. Crowding of the cells and 
nuclei of the epithelial layer may be responsible for this change. As shown in Table 1 , 
one can use the observed changes in tissue scattering to classify tissue quite 
successfully. Light scattering spectroscopy is a procedure that can be used to obtain 
1 0 information about the number and the size of nuclei of the epithelial cell layer. 

Epithelial cell nuclei are the primary targets of reflected light that is singly scattered 
before it is collected by a probe. The intensity and oscillations of this singly- 
backscattered light are characteristic of the number and size of its target nuclei. This 
technique is used to characterize pre-cancerous and early cancerous changes in the 
1 5 colon, the oral cavity, the bladder and Barrett's esophagus (BE). Results of this 
technique are included for the data set of a preferred embodiment to illustrate the 
information that can be acquired and combined with fluorescence and reflectance 
spectroscopies. Light scattering spectroscopy outperforms the other two methods in 
terms of its ability to separate the dysplastic from non-dysplastic BE sites. 
20 The combination of all three techniques provides an extremely sensitive and 

specific tool for the detection of dysplasia in Barrett's esophagus (BE). In this case, 
agreement with histopathology is achieved in terms of separating high-grade dysplasia 
from non-dysplastic and low-grade dysplasia sites. Additionally, all sites are classified 
correctly as dysplastic or non-dysplastic, with the exception of one site. The observed 
25 improvement is expected, since each one of the techniques examines different features 
of tissue biochemistry and morphology that can be altered during the development of 
dysplastic changes. 

Pancreatic carcinoma is one of the first five leading causes of death in Western 
countries and has a very poor prognosis after initial diagnosis. This is due to late 
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presentation of symptoms and the fact that only about 5-20% of all patients are eligible 
for resection. Adenocarcinoma of the pancreas arises from the ductal epithelium. In the 
precancerous stages ductal epithelial cells undergo morphological changes similar to 
those of Barrett's esophagus (BE) including increasing nuclear size and crowding. 
5 Several clinical conditions exist which may allow detection of precancerous or early 
cancerous changes. These include acute or chronic pancreatitis due to focal obstruction 
of the pancreatic duct, acute pancreatitis due to intraductal papillary mucinous tumor of 
the pancreas (IPMT), and patients with a strong family history of pancreatic cancer. 
Detection of dysplasia in this setting is currently based on the detection of stricturing or 
10 dilation of the pancreatic duct, and exfoliative cytology. However, these methods 
cannot reliably distinguish dysplasia and inflammation and have an overall poor 
sensitivity of 44-70%. 

The present invention includes methods for performing tri-modal spectroscopic 
analysis within the main pancreatic duct using the technique of endoscopic retrograde 
1 5 cholangiopancreatography (ERCP), and obtained spectra from the pancreatic duct 
epithelium. 

The majority of patients with ductal adenocarcinoma of the pancreas present 
with late stage tumors that are not amenable to curative therapy. A significant minority 
of patients present for evaluation at a time where early dysplasia of carcinomas are 
20 detectable. These include strictures of the main pancreatic duct (MPD), and IPMT 
presenting as acute pancreatitis, and in patients with a strong family history of 
carcinoma of the pancreas. In all of these cases, the presence of dysplasia, and the 
distribution of dysplasia along the length of the pancreas are critical to management 
decisions (whether to remove part of all of the pancreas). The present invention can be 
25 directed to the use of tri-modal spectroscopy for the detection of dysplasia an invasive 
cancer in the pancreatic duct. The present invention can evaluate ex vivo or in vivo 
specimens of patients to detect pancreatic cancer of dysplasia. Spectra can be collected 
from the MPD at 1 cm intervals within 3 hours of resection and before formalin fixation. 
Tri-modal spectral analysis can be performed to evaluate for components which 
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accurately discriminate histological categories. These spectral algorithms can be 
applied in vivo to patients undergoing endoscopic retrograde cholangiopancreatography 
for evaluation of the 3 conditions mentioned above. Spectra can be collected in vivo at 
1 cm intervals along with entire length of the duct. Spectral signals can be compared to 
5 histology among the patients whose clinical condition dictates surgical removal of the 
pancreas. Anatomic locations of the spectral signals can be matched according to the 
distance from MPD orifice. 

The present invention can be employed for the detection of pre-invasive disease 
of the cervix either alone or at the time of colposcopy following an abnormal Pap smear. 

Q 1 0 The present invention can also be used for a non-invasive method of monitoring the 

JlJ progress of medical therapies for pre-invasive disease. 

W The Pap smear is a screening tool for cervical tissue abnormalities. Abnormal 

fg Pap smears are routinely followed up by colposcopy. This process uses a low-power 

: : i 

^ binocular microscope for the identification of abnormal cervical epithelium, which is 

1 5 subsequently biopsied and examined histopathologically. It is estimated that in expert 

H- hands the sensitivity and specificity of colposcopy are 94% and 48%, respectively. 

S A preferred embodiment of the present invention include spectroscopic 

! y techniques to evaluate cervical epithelium in vivo. Tissue fluorescence spectra excited 

at 337, 380 and 460 nm were acquired during colposcopy from normal and suspicious 
20 sites within the ectocervical and endocervical regions. Suspicious sites were biopsied 
and classified histopathologically. An initial set of spectra was analyzed and statistical 
methods were developed to optimize the agreement between spectroscopic and 
histopathological classification. When these methods were used prospectively to 
classify a second set of data, it was found that squamous intraepithelial lesions (SILs) 
25 can be identified spectroscopically with 82% sensitivity and 68% specificity when 

compared to histopathology. The present invention improves upon the sensitivity and 
specificity of spectral analysis of cervical epithelium in a real-time in vivo approach. 
The method employs tri-modal spectroscopy (TMS), the combined use of intrinsic 
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fluorescence spectroscopy (IFS), diffuse reflectance spectroscopy (DRS), and light 
scattering spectroscopy (LSS). 

The method can compare the spectra obtained using tri-modal processing with 
the histologic diagnosis of the area of epithelium biopsied. This process provides 
5 patterns that are predictive of histologic dysplasia in a prospective fashion, thus 
allowing the clinician to increase the positive histologic lesions of the cervical 
epithelium prospectively. This can be of immense value in following patients on 
clinical trials in order to determine the response to medical treatments of cervical 
dysplasia. 

p 10 Spectra can be acquired from the normal squamous ectocervix and suspicious 

U sites within the transformation zone. The latter may be biopsied immediately following 

j*y spectral acquisition. Data has been collected from 34 patients, 42 normal ectocervical 

i J. 5 

El sites (not biopsied), 15 benign biopsied sites (12 classified as squamous metaplasia and 

3 as mature squamous epithelium) and 10 high-grade squamous intraepithelial lesions 
5 15 (HSILs). 

U Differences in the intrinsic fluorescence intensity and/or lineshape are observed 

q for several excitation wavelengths as shown in Figures 9A-9D. 

1 y Instead of using principal component analysis for the assessment of the 

diagnostic potential of intrinsic fluorescence spectroscopy, the spectra may be 
20 decomposed as a linear combination of the NADH and collagen EEMs extracted from 
the measurements performed during variceal ligation. 

A significant decrease is observed in collagen fluorescence of the benign 
biopsied sites and HSILs compared to that of the normal squamous epithelium as seen 
in Figure 10. These changes result from differences in the levels of expression of 
25 matrix metalloproteinases (MMPs), a class of enzymes responsible for collagen 

degradation. Differences in the levels and/or patterns of expression of MMP-2 have 
been reported between normal squamous epithelium, squamous metaplasia and SILs. 
Additionally, an increase in the NADH fluorescence is noted for the HSILs compared to 
that of benign biopsied sites as seen in Figure 1 1A. This increase can be the result of an 
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increase in the number of epithelial cells and/or their metabolic activity (52). Using 
logistic regression and "leave-one-out" cross-validation, as in the case of Barrett's 
esophagus, HSILs can be distinguished from the benign biopsied sites with a sensitivity 
of 80% and a specificity of 67%. 
5 Analysis of the diffuse reflectance spectra using the periodic component in the 

scattered light indicates that the reduced scattering coefficient of the HSILs is generally 
lower than that of the squamous metaplastic and mature squamous epithelium biopsied 
sites as seen in Figure 1 IB. When the value of ju s ' (A) at 400 run is used to perform 
logistic regression and cross-validation, 7/10 HSILs (70% sensitivity) and 9/1 5 benign 
O 10 biopsied sites (60% specificity) are classified correctly. 

m Extraction and analysis of the light scattering spectra indicated that the variation 

J *i in nuclear size, i.e. the standard deviation of the nuclear size population corresponding 

B to a particular site, has a higher diagnostic value than the percentage of enlarged nuclei 

hi 

s considered in the case of Barrett's esophagus. Using logistic regression and cross 

1 5 validation with the nuclear size standard deviation and the number of nuclei per unit 
M : area as diagnostic parameters, HSILs may be separated from the benign biopsied sites 

O with 90% sensitivity and 67% specificity as seen in Figure 1 IB. 

g Finally, the three techniques, IFS, DRS and LSS, are combined in a preferred 

embodiment in a manner that classifies a particular site according to the diagnosis that is 
20 consistent with at least two of the three methods of analysis. This approach leads to a 

high level of sensitivity and specificity (100% and 80%, respectively) for the detection 

of HSILs from non-HSIL sites biopsied within the transformation zone when compared 

to any one of the techniques alone as seen in Table 2. 
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Table 2 





benign biopsies vs HSIL 




Sensitivity 


Specificity 


Intrinsic Fluorescence Spectroscopy (IFS) 


80% 


67% 


Diffuse Reflectance Spectroscopy (DRS) 


70% 


60% 


Light Scattering Spectroscopy (LSS) 


90% 


67% 


Tri-Modal Spectroscopy (TMS) 


100% 


80% 



Figure 12 illustrates a process sequence 100 of a preferred embodiment of the 
Z invention in which a region of interest within a lumen of a patient is illuminated 102 

and the fluorescent and reflected light is detected 104. The three components 106, 108 
U 5 and 1 10 are used to generate structural and biochemical information 112. This data can 

be analyzed using discriminant analysis 1 14, the components weighted 116, and a 
diagnosis performed 1 18 in real-time. These measurements demonstrate the ability of 
3 spectroscopic techniques to provide useful information for disease classification in a 

non-invasive manner. While each of the techniques can be used for detecting dysplasia 
t 1 0 in Barrett's esophagus, their combination allows the formation of a comprehensive 

U picture of the biochemical and morphological state of tissue. Specifically, 

decomposition of the intrinsic tissue fluorescence EEMs into EEMs of biochemicals 
such as NADH and collagen can provide details about tissue biochemistry. Reflectance 
and light scattering spectroscopy yield morphological information related to the 
15 connective tissue and the epithelial cell nuclei. As this information is free from artifacts 
introduced by tissue excision and processing, it can help advance the understanding of 
the processes that lead to the progression of dysplasia. Software for performing data 
analysis in real-time enables the applicability of these techniques as a guide to biopsy. 
Methods to image regions of interest using this procedure enables large tissue areas to 
20 be studied rapidly. 

An operating environment for the system includes a processing system such as 
data processing system 34, with at least one high speed processing unit and a memory 
system. In accordance with the practices of persons skilled in the art of computer 
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programming, the present invention is described with reference to acts and symbolic 
representations of operations or instructions that are performed by the processing 
system, unless indicated otherwise. Such acts and operations or instructions are 
sometimes referred to as being "computer-executed," or "processing unit executed." 
5 It will be appreciated that the acts and symbolically represented operations or 

instructions include the manipulation of electrical signals by the processing unit. An 
electrical system with data bits causes a resulting transformation or reduction of the 
electrical signal representation, and the maintenance of data bits at memory locations in 
t the memory system to thereby reconfigure or otherwise alter the processing unit's 

10 operation, as well as other processing of signals. The memory locations were data bits 
in are maintained are physical locations that have particular electrical, magnetic, optical, or 

I n organic properties corresponding to the data bits. 

The data bits may also be maintained on a computer readable medium including 
a magnetic disks, optical disks, organic disks, and any other volatile or non- volatile mass 

1 5 storage system readable by the processing using. The computer readable medium 
includes cooperating or interconnected computer readable media, which exist 
□ exclusively on the processing system or is distributed among multiple interconnected 

processing systems that may be local or remote to the processing system. 

Figure 13 is a flowchart of a preferred embodiment of measuring a structure in a 
20 layer of tissue by analyzing data in real-time in accordance with the present invention. 
The real-time analysis may be used to diagnose a pre-cancerous condition in, for 
example, epithelial tissue in the cervix and GI tract. The measurement process begins 
with the step 200 of calibration. The diffuse reflectance spectrum (DRS) is then 
analyzed per step 300. In step 400 the intrinsic fluorescence spectrum (IFS) is analyzed. 
25 The next step in the process includes the step 500 of analyzing the light scattering 

spectrum (LSS). The method then proceeds to step 600, which represents the combined 
analysis for the previous modes of this tri-modal system analysis. 

In a preferred embodiment, during an endoscopic procedure, fluorescence 
spectra at a plurality of wavelengths, for example, at eleven laser excitation wavelengths 
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between 337 and 620 run and one white light (350-750 nm) reflectance spectrum is 
acquired in less than one second. Light delivery and collection is mediated through an 
optical fiber probe. The acquired spectra contain information about the uppermost 
tissue layers, where almost 90% of cancers begin, 
5 From the recorded fluorescence and reflectance spectra, three types of 

spectroscopic information is extracted: IFS, DRS and LSS. The IFS refers to the 
recovery of tissue fluorescence spectra that are free of distortions introduced by tissue 
scattering and absorption. To remove these distortions, the measured fluorescence and 
M reflectance spectra are combined using, for example, a photon-migration-based method. 

□ 10 The extracted IFS is decomposed to provide quantitative information on the 
% \ biochemical tissue composition and the changes that take place in pre-cancerous tissues. 

j£J The measured reflectance spectra consists mainly of photons that are scattered 

hi many times before being detected. A model that is based on diffusion theory is used to 

Pi describe the diffusely reflected light and thus to extract information about the absorption 

1 5 and the reduced scattering coefficients of tissue. The reduced scattering coefficient 
to depends mainly on the morphology of connective tissue, which provides structural 

m support for the epithelium, the most superficial tissue layer. 

The combined use of the IFS, DRS and LSS serves as a tool for biophysical 
tissue characterization and detection of pre-cancerous lesions. In a preferred 
20 embodiment, the TMS analysis may be performed in 2-8 seconds from time of data 
collection. This TMS analysis can provide a real-time guide to biopsies. 

Figure 14 illustrates the details of the step 200 of calibration in accordance with 
a preferred embodiment of the present invention. The process of calibration to remove 
variations in spectra associated with, for example, instruments begins with step 210 of 
25 calibrating any data files in the system that are used to diagnose diseased tissue. The 

next step in the calibration sequence includes the step 220 of calibrating the system with 
a reflectance standard. This standard may include but is not limited to, for example, the 
use of a white light standard such as Spectralon calibration of reflectance. The 
calibration method 200 includes the step 230 of calibrating the system with a 
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fluorescence standard. The fluorescence standard may be, for example, and without 
limitation, Rhodamine. The next step in the calibration method 200 includes the step 
240 of removing the background fluorescence such as, for example, by subtracting any 
background fluorescence. The spectral response of the instrument is calibrated per step 
5 250. 

Figure 15 is a flow chart of a preferred embodiment of the present invention 
detailing the analysis of the DRS. The process includes the generation of a look-up 
table which describes the Mie scattering parameter for different diameters of scattered 
particles per step 310. The creation of the look-up table describing Mie scattering 
O 1 0 reduces processing time for Mie scattering calculations for each iteration. Step 320 
jn includes dividing the reduced scattering coefficient values ju s 1 (X) into a plurality of 

j f: different regions. In a preferred embodiment, there are five regions. This division into 

different regions accounts for the lack of one optimum or universal fit for all ju s f (X) 
values, as there are different values for ju s ' (X) in different regions. Thus, the ju s ' (X) 
15 space is separated into a plurality of regions and a minimization process is performed 
for each region separately. In step 330, the calibrated reflectance spectrum is analyzed 
using a model, for example, the model described in previously referenced U.S. Patent 
No. 6,091,984 and referenced herein after as the Zonios model in each of the different 
regions. The Zonios model accounts for one layer tissue architecture and represents an 
20 iterative process to get a best fit between the results of the model and values from an 
actual recording. A preferred embodiment of the present invention includes accounting 
for a multi-layer tissue architecture such as a two-layer model. In the Zonios model 
light is provided onto a surface of a tissue by a light source. The light reflected back 
from the tissue is a function of wavelength and depends on several parameters such as 
25 the absorption coefficient of tissue, the reduced scattering coefficient and the probe. 
Thus, the reflectance spectrum R(X) = F{ji a (X), ju s ' [l\ r c ) wherein ju a (X) is the 
absorption coefficient of tissue indicative of hemoglobin or betacarotine or more 
generally a biochemical that absorbs light, ju s 1 (X) is the reduced scattering coefficient 
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and r c is the effective light collection area of the probe. A best fit is then selected 
between the model and data per step 340. The best fit is determined by optimizing 
values of parameters, for example, a scaling constant, particle size and density and 
hemoglobin concentration and saturation within a predetermined space. In step 350 a 
5 correction is calculated and applied for high hemoglobin concentrations and/or 
absorption if necessary. Typically the fit is not good when the concentration of 
hemoglobin is high. In a preferred embodiment if the hemoglobin concentration fit 
exceeds a predetermined threshold value, the fit is only performed for a region, for 
example, 450-700 nm, Further, a correction may also be calculated and applied for a 
3 1 0 poor fit as indicated by the shape of the reflectance curve if necessary per step 360. 

Figure 16 is a flow chart detailing the analysis of the IFS in accordance with a 
preferred embodiment of the present invention. It should be noted that biological 
tissues are turbid media in which the predominance of light scattering unavoidably 
entangles the effects of absorption and fluorescence. The interplay of absorption and 
1 5 scattering can substantially distort the measured fluorescence. Therefore, it is important 
to extract the intrinsic fluorescence from the measured tissue fluorescence. The intrinsic 
fluorescence is defined as the fluorescence that is due only to fluorophores, without 
interference from the absorbers and scatterers that are present in the tissue. The greatest 
distortion of fluorescence spectra occurs for emission near 420 nm where hemoglobin 
20 absorbs strongly. Removal of such distortion is particularly important for ultraviolet 
excitation wavelengths, which have been shown to contain important diagnostic 
information. 

In step 410 the input data is analyzed using the information from the reflectance 
fit provided by the Zonios model. The inputs include ju[ (X), ju a (X), fluorescence F(X), 
25 and R obs ; wherein R obs is the observed white light reflectance. Data is also analyzed 
using a method based on photon-migration theory referred to hereinafter as the Zhang 
model. The Zhang model is an improvement over the model described in U.S. Patent 
Number 5,452,723 entitled "Calibrated Spectro graphic Imaging" issued on September 
26, 1995, the entire contents of which being incorporated herein by reference. The 
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15 



20 



Zhang model is described in an article entitled, "Turbidity-free fluorescence 
spectroscopy of biological tissue/' by Zhang et al., and published in Optics Letter, Vol. 
25, No. 19, October 1, 2000, the entire contents of which being incorporated herein by 
reference. The Zhang model is based on the fact that when reflectance and fluorescence 
are collected from a given tissue site with the same probe, absorption and scattering 
distort fluorescence and reflectance spectra similarly, because the diffusely reflected and 
the fluorescence photons traverse similar paths. 

In a preferred embodiment, the fluorescence intensity of the medium can be 
modeled as 



where I x denotes the excitation intensity, <fi xm is the fluorescence quantum yield, in 
which the subscripts x and m denote evaluation of the corresponding quantities at 
excitation frequency v x and emission frequency v m , p ni denotes the fluorescence 



fraction of photons not absorbed during each interaction event and can be considered the 
photon weight reduction factor for that event. Consequently, the expression in brackets 
has the meaning of fluorescence photons per excitation photon. After the photon 
transverses the first i absorption - scattering nodes and experiences propagation loss, 
its weight becomes a l x for each incident excitation photon. For each excitation photon 
arriving at the (i + l)th node, ju fx /(ju^ + ju sx ) photons are absorbed by the fluorophore. 
For each photon absorbed by the fluorophore, (j> xm fluorescence photons are generated. 
The propagation loss of the fluorescence photons is accounted for by the factor a n m l ~ l . 
For a dilute, optically thin (ju f J « l) sample of thickness / , with the same 

amount of fluorophore but without other absorbers and scatterers, one obtains the 
following intrinsic fluorescence, 




(i) 



escape probability for the fluorescence rates, a is 



Ms /(Ma + Ms ) an( * represents the 



f xm =(lJhv x )M f M m hv i 



m 



(2) 
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The fluorescence expression in equation 1 can be expressed in terms of 
reflectance with the help of the following: 

where R is the reflectance of the homogeneous medium of scatterers and absorbers. 
5 The result can be combined with equation 2 to obtain the relationship among intrinsic 
fluorescence spectrum f xm , experimentally measured fluorescence spectrum F xm , and 

reflectance spectrum R m : 

f = $= (4) 



R„ R, 



V ^x^m J 



R 



ox 



R„ 



\ R om 



wherein R 0 is the reflectance in the absence of absorption and g= exp (3 - 1 
10 where /? is a parameter independent of the number of steps n but dependent on the 

scattering properties of the medium, and ju s are the scattering coefficients of the tissue. 

Equation 4 holds for a probe of arbitrary delivery - collection geometry, as long 
as the fluorescence and reflectance are measured with the same probe: note that R Q is a 
function of wavelength. Note that Equation 4 can be used to extract an intrinsic 
1 5 fluorescence excitation - emission matrix (EEM) from an experimentally measured 
fluorescence EEM. When the medium contains multiple fluorophores, one need only 
replace every occurrence of ju^ xm ^ n equations 1 and 2 with a summation of this 
quantity over all fluorophores. Doing so results in the same expression as in equation 4. 
In step 420 the IFS { R obs fit{X)) is obtained. A biochemical decomposition 

20 process may then be performed per step 430. The biochemical decomposition process 
may be a linear decomposition process based on the IFS and provides a quantitative 
measure of the spectra. Further a principal component analysis (PCA) analysis may also 
be performed per step 440. 

Figures 17A and 17B are flow charts detailing the analysis of LSS in accordance 

25 with a preferred embodiment of the present invention. The analysis begins with step 
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510 wherein the observed white light reflectance curve (R obs (X)) and the reflectance fit 

curve [R fit (X)) are smoothed. A moving window average smoothing process may be 

used in a preferred embodiment. An LSS spectrum is generated per step 520 from the 
smoothed data wherein the LSS spectrum is indicative of the difference between the 
5 observed white light reflectance and the reflectance fit 

(LSS spectrum = R obs (X) - R fit (X)) . The LSS spectrum is wavelength dependent. The 

data is then transferred from wavelength (X) space to wavenumber (k) space in step 
530. A first test is then performed per step 540 wherein it is determined if there is a 
O well defined peak in the spectrum in Fourier space. Only if there is a well defined peak 

10a Fourier transform is estimated, otherwise a reliable diagnosis cannot be provided. In a 
preferred embodiment the test for the presence of a well-defined peak may include, for 



i y 



r=maxf R(k) ^ 
^=7^ - 1 pin k 
=min {R(k) ) 



Idk 



If a well-defined peak is present, the process follows to step 550 where a Fourier 
transform (FT) is generated. However, if there is a lack of a well defined peak, the 

15 process returns to step 310. The process then enters step 560 which functions as a 

second test or checkpoint wherein the hemoglobin absorption features are accounted for. 
The second order polynomial is fit to the R obs (X) data. Steps 510, 520 and 530 are then 
repeated if the fit is not accurate. A Fourier transform is then performed to compare the 
resulting FT's from step 550 and step 560. It is then determined if the FT's are the 

20 same. If they are different, the process proceeds to step 570. If the FT's are the same 
however, another hemoglobin or a two layer model, without limitation, is used such as 
described with respect to the DRS analysis in Figure 15. 

In step 570, a third test or checkpoint is performed by analyzing the FT from 
step 550 and ascertaining the presence or absence of a dominant peak. If no dominant 

25 peak is present then the process returns to step 3 1 0 in the DRS analysis as the 

hemoglobin absorption parameter is then probably an artifact. If a dominant peak exists 
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then spectral components other than the dominant peak are removed such as by filtering 
high frequency components per step 574. 

The data is converted to an appropriate scale per step 576 by, for example, 
multiplying by an overall scaling constant to normalize the result. The nuclear size 
5 distribution is determined in step 580. The determination may include the 

determination of the percent enlargement of the nuclear size per step 582. The 
determination of the nuclear size distribution may also include the determination of the 
width of the nuclear size distribution per step 584. Further, determination of the nuclear 

U size distribution may include the determination of the number density of the nuclei per 

S 10 step 586. 

V] Figure 1 8 illustrates a preferred embodiment of a tri-modal system (TMS) 

ill analysis in accordance with the present invention. The TMS analysis includes step 610 

yl of at least one of taking values of diagnostic parameters output from the DRS analysis 

JL (ju a (X) and ju s 1 (X)) , values of diagnostic parameters output from IFS analysis, for 

1 5 example, but not limited to, bio-chemical analysis per step 620 and values of diagnostic 
S3 parameters output from the LSS analysis, for example, without limitation, nuclear size 

5j distribution, size variation, standard deviation and number density per step 630 and 

performing a discrimination analysis in step 640. In a preferred embodiment the 
discrimination analysis may include, for example, but without limitation, Fisher 
20 discrimination analysis which includes assigning weights to each parameter, summing 
the weighted parameters and comparing the sum to a threshold value. The results of the 
TMS analysis are provided as a TMS diagnosis in step 650. The TMS diagnosis 
includes the assessment about the presence or absence of non-dysplastic tissue, high- 
grade dysplastic tissue and dysplastic tissue. 
25 Figure 19 is an alternate preferred embodiment of a TMS analysis in accordance 

with the present invention. In this preferred embodiment a diagnosis from each 
modality of spectroscopy is provided individually and compared before providing an 
overall diagnosis. The values of the diagnostic parameters output from the DRS 
analysis, for example, ///(A), in step 662 are used to provide the diagnosis from the 



{J:\CLIENTS\ip\301505\2020\2020-100\F0201963.DOC;!} 



301505.2020-100 

-31- 

DRS analysis per step 664. For example, the slope and intercept data values are 
determined and compared to threshold values to determine the tissue characteristics. 
Similarly, the values of the diagnostic parameters output from IFS analysis in step 670 
are used to provide a diagnosis commensurate with the IFS analysis in step 672. For 
5 example, the contribution of, but not limited to, NADH and/or collagen, and/or 
porphyrin, to the IFS spectrum are compared to threshold values to provide the 
diagnosis from the IFS analysis. The values of the diagnostic parameters output from 
the LSS analysis in step 680 are in turn used to provide a diagnosis from the LSS 
analysis in step 682. The diagnosis from the LSS analysis includes, for example, the 
h 10 comparison of at least one nuclear size characteristic, such as, nuclear size expressed as 
ri a percent enlargement, number density, size variation, width of distribution versus 

iU crowding to threshold values. Depending upon the tissue and its location, a 

m combination of two or more nuclear size characteristics may be compared to 

corresponding threshold values to provide a diagnosis from the LSS analysis. 
P 1 5 The three diagnoses, specifically the diagnosis from the DRS analysis, the IFS 

k* analysis and the LSS analysis, are compared in step 690. An overall diagnosis is 

2 provided from this comparison in step 692. The overall diagnosis may be derived from 

two similar individual diagnoses, for example, a similar diagnosis from the IFS and LSS 
analyses. Different permutations and combinations of similar diagnoses may be used in 
20 preferred embodiments. 

It will be apparent to those of ordinary skill in the art that methods involved in 
the system of fluorescence, reflectance and light scattering spectroscopy for measuring 
tissue characteristics maybe embodied in a computer program product that includes a 
computer usable medium. For example, such as, a computer usable medium can include 
25 a readable memory device, such as a hard drive device, CD-ROM, a DVD-ROM, or a 
computer diskette, having computer readable program code segments stored thereon. 
The computer readable medium can also include a communications or transmission 
medium, such as, a bus or a communication link, either optical, wired or wireless having 
program code segments carried thereon as digital or analog data signals. 
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While this invention has been particularly shown and described with references 
to preferred embodiments thereof, it will be understood by those skilled in the art that 
various changes in form and details may be made therein without departing from the 
scope of the invention encompassed by the appended claims. 
5 The claims should not be read as limited to the described order or elements 

unless stated to that effect. Therefore, all embodiments that come within the scope and 
spirit of the following claims and equivalents thereto are claimed as the invention. 
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